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Abstract 

Upper bounds at the weak scale are obtained for all XijXim type product couplings of the scalar leptoquark 
model which may affect if" — K" , Bd — Bd and Bs — Bs mixing, as well as leptonic and semileptonic K and 
B decays. Constraints are obtained for both real and imaginary parts of the couplings. We also discuss the 
role of leptoquarks in explaining the anomalously large CP-violating phase in Bs — Bg mixing. 
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The Standard Model (SM), in all probability, is just an effective theory valid up to a scale which is much below 
the Planck scale, and hopefully in the range of a few hundreds of GeV, so that the physics beyond SM can be 
explored at the LHC. Direct production of new particles will definitely signal new physics (NP); while it is an 

■ interesting problem to find out what type of NP is there (commonly known as the 'inverse problem'), it is also 
I well-known that indirect data from low-energy experiments will help to pin down the exact structure of NP, 

including its fiavour sector. The low-energy data, in particular the data coming from the B factories as well as 
from CDF, D0, LHCb (and also from the general purpose ATLAS and CMS experiments) are going to play 
a crucial role in that. There are already some interesting hints; just to name a few [T]: (i) the large mixing 
. , phase in Bg — Bg mixing; (ii) the fraction of longitudinally polarised final states in channels like B — > (t)K* and 

■ B ^ pK*; (iii) the anomalous direct CP-asymmetries in _B ^ ttK decays; (iv) the discrepancy in the extracted 
value of sin(2/3) from Bd J/ipKs and Bd <t>Ks] (v) the larger branching fraction of _B+ t^v compared 
to the SM expectation; and (vi) the discrepancy in the extracted values of Vub from inclusive and exclusive 
modes. While none of them are conclusive proof of any NP, there is a serious tension with the SM when all the 
data are taken together. If all, or most, of them survive the test of time and attain more significance, this will 
indicate a new physics whose flavour sector is definitely of the non-minimal flavour violating (NMFV) type. 

In this paper, as an example of an NMFV new physics, we focus upon the model of scalar leptoquarks (LQ). 
In general, as in any NMFV model, we expect possibly large deviations from the SM in the flavour sector 
observables. LQs that violate both baryon number B and lepton number L must be massive at the level of 
~ 10^^ GeV to avoid proton decay, and are of no interest to us. (There are exceptions; one can construct models 
where LQs violate both B and L and yet do not mediate proton decay. These LQs can be light. For example, 
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see [2].) On the other hand, LQs conserving either B or L or both can be hght, 0(100 GeV), and we discuss the 
phenomenology of only those models that conserve both B and L; one can find extensive discussions on these 
models in [51 m [S]. Vector LQs, as well as some gauge-nonsinglet scalar ones, couple to neutrinos, and their 
couplings should be very tightly constrained from neutrino mass and mixing data. 

Another phenomenological motivation for a LQ model is that this is one of the very few models (R-parity 
violating supersymmetry is another) where the neutral meson mixing diagram gets a new contribution to the 
absorptive part. This, for example, may lead to an enhancement in the width difference AF in the Bg system 
[6] , contrary to what happens in more popular NP models that can only decrease AF [7] . The NP also changes 
the CP-violating phase in Bg — J/ip(j) and hence can help reducing the tension [8] of SM expectation and the 
Tevatron data on the CP- violating phase and width difference for Bg. 

All flavour-changing observables constrain the product of at least two different LQ couplings, one linked with the 
parent flavour and another with the daughter flavour. The product couplings may be complex and it is generally 
impossible to absorb the phase just by a simple redeflnition of the LQ field. We use the data from K'^ — if", 
Bg — Bg and Bd — Bd mixing to constrain the relevant product couplings, generically denoted as AA. For the B 
system, we use the data on AMd,g and the mixing phase sin(2/3d_s), and for the K system, we use AMk and Ek- 
We do not discuss other CP violating parameters like e'/e, since that has large theoretical uncertainties. We 
also discuss the correlated leptonic and semileptonic decays, i.e., the decays mediated by the same LQ couplings. 
While decays to most of the semileptonic channels have been observed, the clean leptonic channels only have 
an upper bound for almost all the cases, except the already observed leptonic decays Kl — >■ e+e", /x+/Lt~. Note 
that the final state must have leptons and hence the bounds are more robust compared to those coming from 
models with only hadrons in the final state. 

A similar exercise have also been undertaken in O |9l [TOl [TTl [12]. We update these bounds with new data 
from the B factories and other collider experiments. In particular, in the subsequent sections, all the previous 
bounds that we quote have been taken from 3 . The — Z?" system has not been considered due to the 
large theoretical uncertainties and dominance of long-distance contributions. We refer the reader to |13j for a 
discussion on the bounds coming from — mixing. Leptonic and semileptonic D and Dg decays have been 
used to put constraints on LQs that couple to the up-type quarks. In particular, LQ contribution might be 
interesting to explain the Dg leptonic decay anomaly [21 [151 IE] ■ The couplings that we constrain are generically 
of the type A^ A*^, where the k-th quark flavour changes to the j-th, but there is no flavour change in the lepton 
sector. One can, in principle, consider flavour changes in the lepton sector too; for an analysis of that type of 
processes, see [14] . However, if one has a vV pair in the flnal state, as in Kl — > tt'^vV, there is a chance that 
lepton flavour is also violated. 

The couplings, which are in general complex, may be constrained from a combined study of CP-conserving and 
CP-violating observables. For neutral mesons, these mean AM as well as ex and sin(2/3d,s)- However, for most 
of the cases, the leptonic and semileptonic decay channels provide the better bound. The analysis has been 
done keeping both the SM and LQ contributions, which keeps the possibility of a destructive intereference, and 
hence larger possible values of the LQ amplitudes, open. 

In Section 2 we briefly state the relevant formulae necessary for the analysis. Section 3 deals with the numerical 
inputs. In Section 4, we take up the analysis, first of the neutral meson mixing, and then the correlated leptonic 
and semileptonic decays. We conclude and summarize in Section 5. 
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2 Relevant Expressions 



2.1 Neutral Meson Mixing 

For the neutral meson system generically denoted by and M*^, the mass difference between the two mass 
eigenstates AM is given by 

-,1/2 



AM = 2Re 



(Mi2 - '-Tr2){M^2 - 1^2) 



AT = -4Im 



(Mi2 - '-T,2){M^2 - ^^12) 



-I 1/2 

(1) 



For the B system, \Mi2\ » |ri2| and AM = 2|Mi2|. For the K system, if the decay is dominantly to the / = 0, 
Imri2 can be neglected and one can write 

_ 1 ImMi2 _ 1 ImMi2 

2V2ReMi2 ^2 AMk ' ^ ' 

Let the SM amplitude be 

Mff = \Mg^ I exp(-2z0sM) , (3) 
where 0sm = for the Bd — Bd system and approximately zero for — and Bg — Bg systems. 

If we have n number of new physics (NP) amplitudes with weak phases 9^, one can write 

n 

M12 = |Mf2^| exp(-2i^sM) + J2 l^i'2l exp(-2i^i) . (4) 

i=l 

This immediately gives the effective mixing phase 9eff as 

. _ 1 , , |Mf,^| sin(2^gM) + E^ \Ml2\ sinim 
- 2 ^'^^'^'^ iMf^^l cos(2^sm) + I^l2l cos(2e,) ' 

and the mass difference between mass eigenstates as 

AM = 2[|Mf2^|2 + ^|Mi'2|'+2|Mf2^|^|Mi2|cos2(^SM-^i) 

i i 

+ 2^^|Mf2||Ml2|cos2(^,-e,)]i/2. (6) 

For the JC" — system, the dominant part of the short-distance SM amplitude is 

^SM ^ mHefflK') ^ p-iy^,v-fm^KfKBKml.S,{x^) , (7) 

where xj = m'j/m'^, fx is the K meson decay constant, and rjK (also called r]cc in the literature) and Bk 
parametrize the short- and the long-distance QCD corrections respectively. The top-quark loop dependent part, 
which is tiny due to the CKM suppression, but responsible for CP violation, has been neglected. The function 
.So is given by 

„ , , 4a; - lla;^ -I- a;^ 3a;^lna; 
^°^")= 4(1 -a;)^ -2(T-r^- 

For the _B° — Bq^ system {q = d for Bd — Bd and q = s for Bg — Bs), we have an analogous equation, dominated 
by the top quark loop: 

^SM ^ = ^{VMrVB^msjlBs^ml^Soix,) . (9) 
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2.2 Leptonic and Semileptonic Decays 



For almost all the cases, the SM leptonic decay widths for neutral mesons are way too small to be taken into 
account, and we can safely saturate the present bound with the NP amplitude alone, except for the Kl sector. 
For example, the branching ratio of Bs is about 3.4 x 10~^ and that of Bd — )■ ^jl^ IJ.~ is about 1.0 x 10~^° 

in the SM, while the experimental limits are at the ballpark of 4-6x10^^. Another exception is the B^ — > l^V 
decay, which proceeds through the annihilation channel in the SM: 

Br(B- ^ ru) = — G^^msmf/llKblVs 1 - ^ , (10) 
where tb is the lifetime of the B meson. 

For the semileptonic decays, we use the following standard convention [17j , given for the B K^*^l+l- 
transition: 

{K{p^)\h,s\B{p,)) = P,F,iq') + (Fo(g2) _ F,{q')) , 

{K*{p2,e)\h^{lTl5>\B{p,)) = ^^q^^I!^e*.q[A^{q')-Ao{q')] 



±ze*(mB + mK-,)A^{q^) T ; P,,(e*. 9)^2(92) 

niB + rriK' 

where P = pi + p2, and q — pi — P2- The pole dominance ensures that ^3(0) = ^o(O), and ^3(9^) can be 
expressed in terms of Ai and ^2- 



2.3 Leptoquarks 

Leptoquarks are colour-triplet bosons that can couple to a quark and a lepton at the same time, and can occur 
in a number of Grand Unified Theories (GUTs) [18], composite models [Ej, and superstring- inspired Eq models 
[9]. In fact, the R-parity violating squarks of supersymmetry, as far as their couplings with fermions go, are 
nothing but LQs. Model-independent constraints on their properties are available ISIIHIS]. 

We focus on the scalar LQ model, which conserves both B and L. The relevant part of the Lagrangian [3] can 
be written as 

-Cs = { (Al5o9l«o-2^l + XRSou'keR)Sl + Aj^e d^eflS'^J + {XlsiUrIl + 
AiJSi 9L«cr2e_R)S'l + Arc, dnlLSl + \LSiql'i'<^2(^°'lL ■ S^'^ } + h.c. 

where {So, So), {Si , Si), and 5" (a = 1, 2, 3) represent the SU{2) singlet, doublet, and triplet LQs respectively. 
A*-' is the coupling strength of a leptoquark to an i-th generation lepton and a j-th generation quark, which is 
in general complex, tr's are the Pauli spin matrices. Note that all the four terms that couple a neutrino with 
a LQ can have potential constraints on neutrino mass and mixing. For example. Si can generate the observed 
neutrino mixing pattern through a type-II seesaw mechanism j^Ul 1211 . However, the constraints also depend 
on the vacuum expectation value of a higher-representation scalar field. That is why we show the non-neutrino 
constraints for these couplings too, keeping in mind that the neutrino constraints may be stronger. 

In this work, we focus only on those processes that involve down-type quarks. Thus, there is no way to constrain 
A_RSo ^-iid Alsi from these processes. In fact, these two sets of coupling can be constrained from processes like 

£)0_£)0 i^jxing and £j — > £j + 7. The latter can be constrained from neutrino mixing too, but as we have just 
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mentioned, the limits would depend on other model parameters. We constrain only five types of scalar LQ 
couplings here: Xlso, ^rSo' '^«-Si , ^ls, ' ^^'^ -^LSi- 

While we have not explicitly shown the generation indices in eq. p^ . it is assumed that the LQs can couple 
with fermions from two di?erent generations. There is another approach which we should mention here. In 
this approach [HIS], one takes the LQ coupling for a single generation, say the third, so that only the third 
generation fermions are a?ected. However, these couplings are taken to be in the weak basis, and when one 
rotates to the mass basis of the quarks, the cross-generational couplings are generated with a mechanism similar 
to that of the Cabibbo-Kobayashi-Maskawa (CKM) mixing. This controls the relative magnitudes, as well as 
the phases, of the LQ couplings, and given a texture in the weak basis, all the couplings in the mass basis are 
present, with predicted magnitudes and phases (the LQ coupling phase also comes from the CKM phase). 

While one gets comparable bounds to that of [3] in this scheme as well, one should note that: 

(i) The mixing matrix for the down-type quarks is not known. What one knows is the misalignment between the 
ul and the bases. Thus, one is forced to consider only the charged-current processes where the misalignment 
(and not the individual rotation matrices) matters. 

(ii) The CKM scheme does not say anything about the rotation matrices for the right-chiral quark sector. 
Whatever one uses there is at best an assumption. A similar analysis has been done for R-parity violating 
supersymmetric models too |22| . 

Thus, we will assume that whatever couplings are nonzero, are so in the physical basis of the quark elds, and 
the phase is arbitrary and not a function of the CKM parameters. 

2.4 Direct Production Limits 

The direct production limits depend on the LQ model, as well as the SM fermions these LQs can couple to. 
The best limits are as follows [23]: rriLQ > 256, 316, 229 GeV for 1st, 2nd, and 3rd generation LQs respectively 
when they are pair produced, and rriLQ > 298, 73 GeV for 1st and 2nd generation LQs when there is single 
production. These are the absolute lower bounds at 95% CL, but the constraints are tighter if the LQ coupling 
to the SM fermions, denoted here by A, is large. For example, if A is of the order of the electromagnetic coupling, 
the first generation LQs were found to have a limit of 275-325 GeV by both the HERA experiments [24]. For 
even larger couplings (A = 0.2-0.5 and above) the LEP experiments exclude a much wider mass range [25], but 
such strong couplings are already almost ruled out if LQ signals are to be observed at the LHC. 

The present-day limits for pair-produced LQs are due to the Tevatron experiments. For a summary, we refer the 
reader to The first generation LQs are searched in 2e-|-2j or le-f 2j-t-MET channel; the second generation 
ones are in 2fi + 2j or lfi + 2j-|-MET channel; and the third generation ones are in 2r + 2b or 26-l-MET channel. 
We refer to the original papers [27] for the details of the analysis. 

The production of LQ states, either single, associated with a lepton (from qg — >■ LQ + i), or in pair, from 
QQi 99 ~^ LQ + LQ, has been studied in detail; for example, the reader may look at [HJ [551 [521 ISO]- At y/s = 14 
TeV, the cross-section of pair production of scalar leptoquarks is about 1 pb [30]. While this goes down 
significantly for the initial LHC run of -y/s — 7 TeV, one expects to see LQ signals upto ttilq = 500 GeV even 
with 5 fb^^ of luminosity. The cross-section for single production depends on the value of A. For A = v47ra, the 
electric charge, the cross-section for LQ plus charged lepton production is about 100 fb for itilq ~ 500 GeV. 
The cross-section is proportional to A^, so we expect events even with ^/s = 7 TeV and A ~ 0.05. Obviously, 
for smaller values of A, pair production will be more favoured, and we expect the preliminary run of LHC to 
establish a limit of the order of 500 GeV. 

In this analysis, we will use a somewhat conservative reference mass value of 300 GeV for every LQ state, 
independent of the quantum numbers and generation. The bounds on the product couplings scale as mj^Q, so 
the bounds that we show should be multiplied by (300/miQ)^. 
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Figure 1: Leptoquark contributions to Af" — Af" mixing. 



2.5 Constraints from Meson Mixing 

Consider the neutral meson M'^ = qjOk- The oscillation can have a new LQ mediated amplitude, with i-type 
leptons and some scalar LQs in the box, as shown in Fig. ((!]). The amplitude is proportional to (A*^,Ay)^. We 
consider, as in the standard practice, a hierarchical couphng scheme, so that we may consider only two LQ 
couphngs to be nonzero at the most. Also, we consider any one type of LQ to be present at the same time. 
This keeps the discussion simple and the numerical results easily tractable; however, this may not be the case 
where we have some high-energy texture of the couplings and there can be a number of nonzero couplings at 
the weak scale. 

For the LQ box, one must consider the same type of lepton flowing inside the box if we wish to restrict the 
number of LQ couplings to 2. The effective Hamiltonian contains the operator Oi, defined as 

dl=[h''PRd]^\bj^PRd]^ , (13) 

(where the subscript 1 indicates the SU(3)c singlet nature of the current), and is given by 

(\*feAig)^ 



'H-LQ 



1287r2 




Oi , (14) 



where ci = 1, C2 = for So,Sq and Si,ci = C2 = l for Si, and ci = 4, C2 = 1 for Si. Therefore, if we are 
allowed to neglect the SM, the limits on the product couplings for {XlSoj ^rSq' ■^RSi/2)j ^ls-i/^' ^'^'^ ^LSi should 
be at the ratio of 1 • ■ ■ The operator Oi is multiplicatively renormalized and the LQ couplings are those 
obtained at the weak scale. The function /(x), defined as 



l-a;2 + 2a;logx 



is always very close to /(O) — 1; note that we have taken all LQs to be degenerate at 300 GeV. 



2.6 Constraints from Leptonic and Semileptonic Decays 

The LQ couplings which may contribute to — , Bd — Bd and Bs — Bs mixing should also affect various 
LQ-mediated semileptonic {h d{s)l'^l~, s — )■ dl^l~) and purely leptonic (-B^j^-) — > , l^l^) decays. 

The estimated BRs of leptonic flavour conserving AB{S) = 1 processes within SM are very small compared to 
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Interaction 


4-ferniion vertex 


Fierz-transformed vertex 




G{dLVL){vLdl) 


hG{dLl^dl){vLl^,VL) 






\G{d''j^Yd''R){eRl,.eR) 




G{dLeR){eRdL) 


■^G(dLJij.dL){eR'y^eR) 


= ^ 


GidRVLWhdR) 
G(dReL){eLdR) 


^G(dRrdR){VLi^J^L) 
^G{dR-f''dR){eLjf,eL) 




G{VLdl) 
2G{dleL){eLdl) 


\G{dll^^dl){VLl,VL) 

Gidlrdl){eLlf.eL) 



Table 1: Effective four-fermion operators for scalar leptoquarks. G generically stands for A^/m| 



their experimental numbers or upper bounds, except for — > e^e^,/^^/i^. Therefore it is quite reasonable 
to ignore the SM effects for these channels while constraining the LQ couplings. For these mixing correlated 
decays, the final state leptons must be of the same flavour. The leptonic decay modes are theoretically clean 
and free from any hadronic uncertainties. The semileptonic modes have the usual form-factor uncertainties, 
and the SM contribution cannot be neglected here. 

To construct four-fermion operators from A type couplings which mediate leptonic and semileptonic B and K 
decays, one needs to integrate out the LQ field. The effective 4-fermi Hamiltonians and vertices which are 
related to the mixing is given in Table [TJ The vertices show that the limits coming from leptonic or semileptonic 
decays will be highly correlated. For charged leptons in the final state, one can constrain RSq, RSi, LSi, or 
LSi type LQs. The bounds for the first three will be the same, which is just twice that of LSi. Similarly, if we 
have neutrinos in the final state, LSq, LSi, or LSi type LQ couplings are bounded, all limits being the same. 

The product LQ coupling may in general be complex. If we neglect the SM, there is no scope of CP violation and 
the data constrains only the magnitude of the product, so we can, if we wish, take the product to be real. In fact, 
if we assume CP invariance, Ks decay channels to constrain only the real part of the product couplings, 

and Kl constrains the imaginary part. For the processes where the SM contribution cannot be neglected, 
we have saturated the difference between the highest experimental prediction and the lowest SM expectation 
with an incoherently summed LQ amplitude. For a quick reference, the data on the leptonic and semileptonic 
channels is shown in Tables [5] and [H Note that these bounds are almost free from QCD uncertainties except 
for the decay constants of the mesons, and hence are quite robust. There are other semileptonic channels which 
we do not show here, e.g., B K*vT', because they yield less severe bounds. 



Mode 


Branching ratio 


Mode 


Branching ratio 


Ks e+e~ 
Kl e+e" 
Bd fJ-'^^J-^ 
Bd ^ T+T- 

Bg /Lt+/X" 


< 9 X lO-'J 
9tl X 10-12 

< 1.0 X 10-^ 

< 4.1 X 10-3 

< 3.3 X 10-8 


Ks fi+fi- 
Kl fJ-'^f-i^ 
Bd M^/U" 
Bs e+e- 


< 3.2 X 10-^ 
(6.84 ±0.11) X 10-^ 

< 1.0 X 10-8 

< 5.4 X lO-'"^ 



Tabic 2: Branching ratios for some leptonic decays of K and B mesons [23]. The limits are at 90% confidence 
level. The SM expectation is negligible. 

The constraints coming from the decay M^{= qjQi^) — > iili can be expressed as 

\X,,X*,\<2^mlQ (16) 
for RSa, RSi, and LSi types, and without the factor of 2 on the righthand side for the LSi type LQs. Here 

Fm = -^Br {M° ^ , 
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Mode 


Branching ratio 


SM expectation 


Mode 


Branching ratio 


SM expectation 






S.Oli:^ X lO-'' 


2.1 X 10"'" 






2.9li;^ X 10-"* 


4.8 X 10""' 


Kt 


— 7r^e"^e^ 


< 2.8 X 10"^° 


2.4 X 10"" 


Kt — 


y TV^ fl'^ fl^ 


< 3.8 X 10-"' 


4.4 X 10-^^ 


Kt 


— y vly 


< 6.7 X 10"** 


(2.8 ±0.6) X 10"" 


K+ 


- — y TT^ 1^17 


(17.3liJj) X 10-" 


(8.5 ±0.7) X 10"" 


K+ - 


y 7r~^6~^C 


(2.88 ±0.13) X 10-^ 


(2.74 ±0.23) X 10"^ 


K+ - 




(8.1 ± 1.4) X 10-* 


(6.8 ±0.6) X 10"* 




7r"e+e" 


< 1.4 X 10"' 


3.3 X 10"" 


^d 




< 1.8 X 10"' 


3.3 X 10"° 


Bd - 


■> K°e+e- 


{l.itW) X 10"^ 


2.6 X 10"'^ 


Bd 




X 10"^ 


(3.3 ±0.7) X lO"'^ 


Bd- 


> K*tJ,+ fi- 


(1.06 ±0.17) X 10"^ 


(1.0 ±0.4) X 10"'' 


Bd ^ 


K*e+e- 


1.39 X 10"® 


(1.3 ±0.4) X 10"® 


Bd 


— 5- tt'^uV 


< 2.2 X 10-* 


(8.5 ±3.5) X 10"* 


Bd - 


K°vV 


< 1.6 X 10"'' 


(1.35 ±0.35) X 10"^ 


Bd 


K*%V 


< 1.2 X 10"* 


3.8 X 10-'' 


B+- 




< 8.0 X 10-® 


(2.03 ±0.23) X 10"* 


B+- 


Tr^/i^/i^ 


< 6.9 X 10"® 


(2.03 ±0.23) X 10"** 


B+ 


— TY^UV 


< 100 X 10"® 


(9.7 ±2.1) X 10"® 


B+ - 


-> K+e+e- 


< 1.25 X 10-^ 


6.0 X 10"'^ 


B+ 




< 8.3 X 10-® 


6.0 X lO"'^ 


B+ 


^ K+'-uV 


< 8.0 X 10"^ 


(12.0 ±4.4) X 10"® 


B+ - 


K+vV 


< 14 X 10"® 


(4.5 ±0.7) X 10"® 


Bs - 




(1.44 ±0.57) X 10"*^ 


1.6 X 10"" 


Bs 


— >■ (j)uT' 


< 5.4 X 10"^ 


(13.9 ±5.0) X 10"® 



Table 3: Branching ratios for some semileptonic K and B decays [231 [211 [311 13S1 [M]- The hmits are at 90% 
confidence leveL Also shown are the central values for the SM. For the SM expectations shown with an error 
margin, we have taken the lowest possible values, so that the LQ bounds are most conservative. The systematic 
and statistical errors have been added in quadrature. 



Gm = /mo tmo MIjo me^l- 4-^1 ^ (17) 

T and /jvfo being the lifetime and the decay constant of M° respectively. Note that Kl has a lifetime two orders 
of magnitude larger than that of A'5 and hence the bounds coming from Kl decays are going to be tighter by 
that amount. 



3 Numerical Inputs 

The numerical inputs have been taken from various sources and listed in Table [H We use the BSW form factors 
|17j with a simple pole dominance, and the relevant form factors at zero momentum transfer = are taken 
as follows [55] : 

Fo^^^(O) = ^^'(0) = 0.38 , ^^o^^"(0) = ^"(0) = 0.33 , 

K* : V{0) = 0.37 , Ai(0) = ^3(0) = 0.33 , Ao(0) = 0.32 , (18) 

while we take ^^^^^''(O) = 0.992. This is not incompatible with the lattice QCD resuh of 0.9560(84) [36]. The 
theoretical uncertainty comes mostly from the form factors, but is never more than 10% for the LQ coupling 
bounds. The bounds are not a sensitive function of the exact values of the form factors, and remain more or 
less the same even when one uses the light-cone form factors. 

The mass differences AM are all pretty well-measured; for consistency, we use the UTfit values [37] ■ We use 
sin(2/3ci) as measured in the charmonium channel [T]. The SM prediction is taken from the measurement of the 
UT sides only since that is least likely to be affected by new physics. (However, this need not be true always. 
For example, if there is a new physics contributing in the Bd — Bd mixing amplitude, the extracted value of Vtd 
may not be equal to its SM value.) For f3s, which is defined as arg(— Vfs V^J/VcsV'^*},), the errors are asymmetric: 

= (o.47t°i?) u {imtllD , (19) 

which we show in a symmetrized manner. The decay constants Jb^ , are taken from |38| as a lattice average of 
various groups. The same holds for JbVBb and ^, defined as ^ = /b^ ■\/ Bb^ / Jb^ \J Bb^ , whose value we take 
to be 1.258 ±0.020 ±0.043. 
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Observable 


Value 


Observable 


Value 


AMk 
AMb, 
sin(2/3d)exp 


5.301 X 10"^ ps^^ 
(0.507 ± 0.005) ps-i 
(17.77 ±0.12) ps-i 
0.55 ±0.01 
0.668 ±0.028 
0.731 ± 0.038 
(0.43 ±0.17) U (1.13 ±0.17) 


\eK\ 
Bk 

VBk 
fK 

fB. 
JBs I fBi 
IBs \J Bbs 


(2.228 ±0.011) X 10"^ 
0.75 ±0.07 
1.38 ±0.53 
160 MeV 
(228 ± 17) MeV 
(1.199 ±0.008 ± 0.023) 
(257±6±21) McV 



Table 4: Input parameters. For the form factors, see text. 



4 Analysis 

4.1 Neutral Meson Mixing 

While our bounds are shown in Table [5] following the procedure outlined in Section [2751 let us try to understand 
the origin of these bounds. 

Take Figure [2] (a) as an example, which shows the bounds on the real and imaginary parts of AiiA*2. This is 
shown for the triplet LQ 5*1; all LQs produce a similar diagram, with the limits properly scaled. To get an idea 
of the scaling, one may again look at Table EJ and scale accordingly. 

For the K system, we use AMk and {ekI as the constraints. The SM part is assumed to be dominated by 
the short-distance contributions only. Note the spoke-like structure; this is because gives a very tight 
constraint on Im(Mi2) and only those points are chosen for which (AA*)^ is almost real. However, as we will 
see later, all the bounds except those for the LSq type LQs will be superseded by those coming from leptonic 
and semileptonic K decays; however, j = 3 bounds will stand. 

A similar analysis is shown for the Bd — Bd system in Figure [2] (b) and Table \E\ Note that the bounds on the 
real and the imaginary parts of any product coupling are almost the same. This is, of course, no numerical 
accident. To understand this, let us analyse the origin of these bounds. There are two main constraints for the 
Bd system: AMd and sin(2/3d). There will be a region, centred around the origin of i?e(AA) — /m(AA) plane 
(since AMd can be explained by the SM alone), where |AA| is small and the phase can be arbitrary. At the Icr 
level, this region appears to be small, because the measured value of sin(2/3d) from the charmonium channels 
is just barely compatible with that obtained from a measurement of the sides of the unitarity triangles. The 
region expands if we take the error bars to be larger. This is the SM-dominated region, where LQ creeps in 
to whatever place is left available. Any analysis, taking both SM and LQ but assuming incoherent sum of 
amplitudes, should generate this region only. 

However, there is always scope for fully constructive or destructive interference between SM and any NP. 
Consider a situation where the LQ contribution is large, so large that even after a destructive interference with 
the SM amplitude, enough is left to saturate AMd- This LQ-dominated region (this is true for all NP models 
in general) gives us the bounds, and in the limit where the SM can be neglected, the bounds on i?e(AA) are 
almost the same as on /m(AA). The alignment of the fourfold symmetric structure is different from Figured] 
(b) because of the sizable value of sin(2/3d). 

The limits for the Bg system are shown in Figure [3] Note that the origin is excluded at the la level; this is due 
to the large observed values of (3s: Ps = (25 ± 10)° U (65 ± 10)° in the first quadrant and a mirror image in the 
second quadrant. 

The magnitude of the product is bounded to be less than 0.08 at the Itr level for the triplet LQ, and scaled 
according to Table[5j For i = 2, the relevant coupling mediates the leptonic decay Bg — > and semileptonic 
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Process 
& indices 


LQ Type 


Real 
Only 


Real part 
of Complex 


Img part 
of Complex 


|AA*| 


(zl)(z2)* 


LSq, RSq, RSi 

LSi 
Lsl 


0.008 
0.0055 
0.0036 


0.008 
0.0055 
0.0036 


0.008 
0.0055 
0.0036 


0.008 
0.0055 
0.0036 


Bd — Bd 


LSq, RSq, RSi 

LSi 
Lsl 


0.009 
0.0063 
0.004 


0.022 
0.016 
0.010 


0.022 
0.016 
0.010 


0.027 
0.019 
0.012 


Bs — Bs 


LSq, RSq, RSi 

LSi 
LS\ 


0.05 
0.034 
0.02 


0.13 
0.09 
0.06 


0.13 
0.09 
0.06 


0.18 
0.13 
0.08 



Table 5: Bounds from the neutral meson mixing. The third column shows the bounds when the couplings are 
assumed to be real. The last three columns are for complex couplings. 




Figure 2: (a) Allowed parameter space for AiiA*2 for \lSi type couplings, (b) The same for AiiA*3. 
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Kl{S) Decay 


Coupling 


|AA*| 


Bd(s) Decay 


Coupling 


|AA*| 


Ks e+e' 
Ks M^M" 
Kl -> e+e~ 
Kl -J> 


(12)(11)* 
(22)(21)* 
(12)(11)* 
(22)(21)* 


1.8 X 10-^ 
5.5 X 10^3 
2.4 X 10-"* 
6.4 X 10-6 


Bd T+T- 
Bs tJ.^tJL'' 


(21) (23)* 
(31)(33)* 

(22) (23)* 


2.8 X 10-^ 

1.2 X IQ-i 

4.3 X 10-3 



Table 6: Bounds from the correlated leptonic Kl(S) and Bd(s) decays. The LQs are either of RSq, RSi, or LSi 
type. For LSi type LQ, the bounds are half of that shown here. 



Bd — i4r(*)/i+/^- decays. We will see in the next part, just like the K system, that the constraints coming from 
such decays are much stronger. The same observation is true for i = 1. Again, only for XlSq type couplings, 
there is no leptonic or semileptonic contributions (the down-type quark current couples with the neutrino current 
only), and the bounds coming from the mixing stand. Thus, for j = 1,2 and any other LQ except 5*0, it is 
extremely improbable that the LQ contribution explains the large mixing phase. 

What happens for i = 3? This will mediate the decays Bg — > r+r" and B XgT+r' . While there is no data 
on these channels yet, we may have a consistency check with the lifetime of Bg- This tells us that couplings 
as large as 0.05 are allowed, but the decay Bg — >• t+t~ should be close to the discovery limit. This will be an 
interesting channel to explore at the LHC. There is an exception: if we consider XlSo type couplings, neutrinos 
flow inside the box, and then we have final-state neutrinos, and not t leptons. 

Note that the box diagram with leptoquarks and leptons has a nonzero absorptive part, which is responsible for 
the corresponding correlated decays. This affects the width differences /^Td,s- As has been shown in [6], NP that 
contributes to AF may enhance the mixing phase in the B^ — Bg box, contrary to the Grossman theorem [7], 
which tells that the mixing phase in the Bg system must decrease due to NP if there is no absorptive amplitude 
in the box diagram. The effect on AFd/Fd is negligible; with the bounds that we get here, it is never more than 
1%, or even less (note that [11 uses a LQ mass of 100 GeV and we need to scale their results). For Bg, AFg/F,, 
may go up to 30% without significantly enhancing the leptonic branching ratios like Bg — > r+r^, and one can 
also get a significant nonzero phase in the Bg — Bg mixing that is indicated by the present experiments [8]. 




(a) (b) 

Figure 3: (a) Allowed parameter space for Ai2A*3 (b) The reach for the angle fig. For more details, see text. 
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Decciy clia/iinGl 


Coupling 


|AA*| 


Decay channel 


Coupling 


|AA*| 


i\g — y IT c e 


(11) 


(12)* 


2.8 X 10" 


-'6 




(llj(12j 


2.8 X 10" 


-b 


i\S — r IT ^ IX 


(21) 


(22)* 


4.6 X 10" 


-3 


Kl - 


+ — 

> 7T fl^ 


\Zl )[ZZ ) 


5.9 X 10" 


-5 


/ V — T TT 6 e 


(11) 


(12)* 


1.2 X 10" 


-3 


K+ — 


1 1 


('91 Voo^* 

\Z1 )\ZZ j 


9.5 X 10" 


-4 


_ V. /rrO 7 /T7 

J \ ~r TV Jyiy 


(*1) 


{i2)* 


4.6 X 10" 


-4 












— ?■ TT e e 


(11) 


(13)* 


1.1 X 10" 


-3 


T-»-l- 

B+ - 


1 1 


(11)(13) 


5.0 X 10" 


-4 


iid — >■ TT /U^/i 


(21) 


(23)* 


1.2 X 10" 


-3 


n+ 

13^ — 




(zi j(zo j 


4.6 X 10" 


-4 


±>(i — )■ A e e 


(12) 


(13)* 


3.6 X 10" 


-4 


B+ ^ 


K+e+e- 


(12)(13) 


7.0 X 10" 


_3 




(12) 


(13)* 


9.7 X 10" 


-4 


Bd^ 






1.1 X 10" 


-3 




(22) 


(23)* 


1.5 X 10" 


-3 


B+ ^ 




(22)(23)* 


5.6 X 10" 


-3 


— > 7r''z^77 


(^l) 


{ay 


4.4 X 10" 


-2 


B+ 


— )• TT^iyV 


(*1)(*3)* 


2.0 X 10" 


-2 


Bd K"uV 


{i2) 


[ay 


2.6 X 10" 


-2 


B+ - 


K+vV 


(i2)(i3)* 


6.5 X 10" 


-3 


Bd K^*vV 


(i2) 


[ay 


1.5 X 10" 


-2 


B+ - 


^ K+*vV 


(i2)(i3)* 


1.2 X 10" 


-2 


Bs 1-'''' 


(22) 


(23)* 


7.9 X 10" 


-4 


Bs 


— > (jyvV 


(i2)(i3)* 


9.1 X 10" 


-2 



Table 7: Bounds from the correlated semileptonic B and K decays. The LQs are either of RSq, RSi, or LSi 
type. For LSi type LQ, the bounds are half of that shown here. For the final state neutrino channels, the LQ 
can be LSq, LSi, or LSi type, all giving the same bound. 



4.2 Leptonic and Semileptonic Decays 

We have assumed only two LQ couplings to be present simultaneously, with identical lepton indices. Thus we 
will be interested only in lepton flavour conserving processes. A similar analysis was done in j39, for vector 
LQs. Our bounds are shown in Table [Bland Table [71 

Apart from the leptonic decays, the SM amplitudes can be neglected as a first approximation. Thus, one 
may saturate the experimental bounds with the LQ amplitude alone. This generates most of the numbers in 
Table [6] For Kl decays, we consider the SM part too, and add the amplitudes incoherently. Note that 
decays only constrain the imaginary part of the LQ coupling. This can be understood as follows. Consider the 
decay Kl — > /^^A* • While in the limit of CP invariance, one can write Kl = [K^ — K'^)/^/2, it is A21A22 that 
mediates K'^ decay and A21A22 that mediates K^ decay. Taking the combination, the imaginary part of the 
coupling is responsible for Kl decays, and the real part is responsible for Ks decays. As mentioned earlier, 
Bs T^T^ does not have a limit yet, but the SM expectation is about 0(10^^), and if IA32A33I ^ 10^^, one 
expects the BR to be of the order of 4 x 10""^. 

For the channel K^ — > tt^vv, the outgoing neutrino can have any flavour, and so the bound is valid for i = 1, 2, 3. 
However, these bounds are valid when one can have a neutrino in the final state, i.e., for LQs of the L category, 
which couple with lepton doublets. 

Semileptonic decays give the best bounds, but they are the least robust one, considering the uncertainty in the 
form factors. While we take the BSW form factors [17], the lattice QCD or light-cone sum rules based form 
factors may change the final results by at most 10%. To be conservative, we saturate the difference between 
the SM prediction and the maximum of the data by LQ contributions. 

Let us just say a few words about B^ — > t^V. In the SM, the branching ratio can be worked out from eq. 
(fTO]) and is (9. 312:3) x 10"^ where the major sources of uncertainty are \Vub\ and fs- The observed number, 
(14.3 ± 3.7) X 10~^ [T] is a bit above the SM prediction. The tension can be alleviated with LSq or LSi type 
leptoquarks; the necessary combination is A31A33, and the bounds that we have obtained on this particular 
combination in Table [9] can easily jack up the branching ratio to the observed level. A similar exercise has been 
done for the leptonic Dg decays in [TSj. 

We have summarized our bounds in Tables [51 [H and [TOl These tables contain no new information, but just 
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LQ 
type 


indices 


Previous 
Bound 


This analysis 


From Mixing 


From Decay 


Real part 


Imag. part 


Channel 


Bound 




(*l)(z2)* 


1.8 X 10-'* 


8 X 10-3 


8 X 10-3 


Kl — > 7r'V77 


4.6 X 10-4 


RSq, 

RSi/2 


(11)(12)* 


2.7 X 10-3 


8 X 10-3 


8 X 10-3 


K+ TT+e+e- 
Kl 7r°e+e- 


1.2 X lO-'' 
{2.8 X 10-^) 


(21)(22)* 


5.4 X 10"^ 


8 X 10-3 


8 X 10-3 


K+ -> TT+fl+fl- 

Kl H+fi- 


9.5 X 10-4 
6.4 X 10-^ 


(31)(32)* 


0.018 


8 X 10-'' 


8 X 10-'' 


— 


— 


LSi/2 


(11)(12)* 


1.8 X 10"-* 


5.6 X 10-3 


5.6 X 10-3 


K^ — >• TT+e+e- 
Kl 7r°e+e- 


1.2 X iO-^ 
{2.8 X 10-^) 


(21)(22)* 


1.8 X 10"^ 


5.6 X 10-3 


5.6 X 10-3 


K^ — > 7r+/i+^- 
Kl M+A*- 


9.5 X 10-^ 
{6.4 X 10-^) 


(31)(32)* 


5.4 X 10"^ 


5.6 X 10-3 


5.6 X 10-3 


Kl -> TT^h'T' 


4.6 X 10-4 


LSi 


(11)(12)* 


1.8 X 10-'* 


3.6 X 10-3 


3.6 X 10-3 


K+ — > TT+e+g- 
Kl 7r°e+e- 


6.0 X 10-4 
{1.4 X 10-^) 


(21)(22)* 


2.7 X 10"'^ 


3.6 X 10-3 


3.6 X 10-3 


K+ 7r+/x+^- 

Kl fJ-'^^J-- 


4.8 X 10-4 
{3.2 X 10-^) 


(31)(32)* 


1.8 X 10-" 


3.6 X 10-3 


3.6 X 10-3 


Kl — > 7r'V77 


4.6 X 10-4 



Table 8: Bounds coming from K^ — K'^ mixing and correlated decays. The better bounds have been emphasized. 
Note that Ks decays constrain Re (AiiA*2) while K^ decays constrain only the magnitudes; however, in view 
of a tight constraint on the imaginary part, the bound from K^ decay can be taken to be on the real part of 
the product coupling. Here and in the next two tables, all numbers in the "Previous bound" column are taken 
from [3], with scaling the LQ mass to 300 GeV. 



LQ 

type 


indices 


Previous 
Bound 


This analysis 


From Mixing 


From Decay 


Real part 


Imag. part 


Channel 


Bound 




(*1)(*3)* 


0.036 


0.022 


0.022 




2.0 X 10-''^ 


RSq, 

RSi/2 


(11)(13)* 


0.054 


0.022 


0.022 


B+ TT+e+e- 


5.0 X 10-4 


(21)(23)* 


7.2 X 10-3 


0.022 


0.022 


B+ TT+n+n- 


4.6 X 10-4 


(31)(33)* 


0.054 


0.022 


0.022 


Bd -> T+T- 


1.2 X 10-* 


LSi/2 


(11)(13)* 


0.054 


0.016 


0.016 


B+ — > TT+e+e- 


5.0 X 10-4 


(21)(23)* 


7.2 X 10-3 


0.016 


0.016 


£?+ 7r+^+/x- 


4.6 X 10-4 


(31)(33)* 


0.054 


0.016 


0.016 


_B+ — > TT^iyp 


2.0 X 10-^ 


LSi 


(11)(13)* 


0.036 


0.010 


0.010 


B+ TT+e+e- 


2.5 X 10-4 


(21)(23)* 


3.6 X 10-3 


0.010 


0.010 


B+ TT+n+n- 


2.3 X 10-4 


(31)(33)* 


0.027 


0.010 


0.010 


Bd T+T- 


6.2 X 10-'^ 



Table 9: Bounds coming from Bd — Bd mixing and correlated decays. The better bounds have been emphasized. 
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LQ 
type 


indices 


Previous 
Bound 


This analysis 


Prom Mixing 


From Decay 


Real part 


Imag. part 


Channel 


Bound 


LSo 


(i2)(i'S)* 


0.36 


0.13 


0.13 




6.5 X 10^'^ 


RSl/2 




5.4 X lO^'' 


n 1 

u. ±o 




Rj — i R'Op+p- 


9 /i V m^'i 

O . U A lU 




7.2 X 10^'^ 


0.13 


0.13 


R . _i. K*n~^n~ 


1.1 X 10~''' 




.09 


0.13 


0.13 






LSi/2 




5.4 X lO""* 


0.09 


0.09 


Bd K^e+e- 


3.6 X 10-"^ 


(22)(23)* 


7.2 X 10--'' 


0.09 


0.09 


Bd ^ K*fi+fi- 


1.1 X 10-^ 


(32)(33)* 


0.054 


0.09 


0.09 


B+ K+vV 


9.3 X 10-^ 


LSi 


(12)(13)* 


2.7 X 10--' 


0.06 


0.06 


Bd ^^"6+6- 


1.8 X 10-'' 


(22)(23)* 


3.6 X 10"^ 


0.06 


0.06 


Bd ^ K*ii+ix- 


5.5 X 10-^ 


(32)(33)* 


0.045 


0.06 


0.06 


B+ K+uV 


6.5 X iO"^ 



Table 10: Bounds coming from Bg — B^ mixing and correlated decays. The better bounds have been emphasized. 

shows the best bound for a given LQ type and a given set of indices. They further show that 

(i) Except for R-type LQs with indices (31) (33), semileptonic, and in a few cases leptonic, decays give the best 
constraints. In most of the cases they are one or more orders of magnitude stronger than those coming from 
the mixing, so with those LQs, one should not expect much discernible effects from CP asymmetries. 

(ii) While the bounds coming from decays are only on the magnitude of the product couplings, information on 
the complex weak phases of these couplings must come from mixing data, unless one makes a careful study of 
semileptonic CP asymmetries. 



5 Summary and Conclusions 

In this paper we have computed the bounds on several scalar leptoquark coupling combinations coming from 
M° — M° mixing as well as leptonic and semileptonic decays. Though such an analysis is not new, we have 
implemented several features in this analysis which were not been taken into account in earlier studies. Apart 
from the improved data on the B system, we have also used the data on CP violating phases, and obtained 
nontrivial constraints on the real and imaginary parts of the couplings. We note that for the gauge-singlet LQ 
with XlSo type couplings, it is possible to alleviate the mild tension between the measured and predicted 
values of sm{2/3d), as well as to explain the large mixing phase in the Bg system. For this type of LQs, there 
are no modifications in leptonic or semileptonic channels, unless we consider final-state neutrinos. 

For all other type of LQs, leptonic and semileptonic decays provide the better constraints (the exceptions are 
final-state t channels). We do not expect any effects on nonlcptonic final states like those coming from, say, 
R-parity violating supersymmetry with A' type couplings. While the bounds coming from the leptonic channels 
are quite robust (apart from the mild uncertainty in the meson decay constants), those coming from semileptonic 
decays have an inherent uncertainty of the order of 10-15%, whose origin is the imprecise nature of the form 
factors. 
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